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Q>^ ■ A review is presented of what we understand of colour confinement in QCD. Lattice 

formulation provides evidence that QCD vacuum is a dual superconductor: the chromo- 

electric field of a qq pair is constrained by dual Meissner effect into a dual Abrikosov flux 

>• I tube and the static potential energy is proportional to the distance. 
O 

O ' 1- Introduction 
l> . 

O ■ Lattice formulation is a gauge invariant regulator of non abelian gauge theories. 
• Numerical simulations on the lattice produce from first principles regulated correlators of 
physical quantities. 

Simulations can be used to compute quantities involving low energy modes, which are 
out of reach of perturbation theory, such as weak interaction matrix elements, masses, ma- 
trix elements of operators in the light cone expansion. The typical problems encountered 
in this "phenomenological" use of lattice are the removal of the cut off (renormalization) , 
and limitations in computer power. 
^ ! Simulations can also be used to test theoretical ideas and to investigate the structure of 
I the theory. An example of investigation of "theoretical" type is the study of the mechanism 
of confinement. The typical difficulty in this approach is to have good theoretical ideas 
to test numerically, and possibly to falsify from the first principles. 

Apart from confinement itself there are a few fundamental issues at the background of 
our understanding of QCD. Among them 

1) The limit. The conjecture^ is that the basic properties of a gauge 

theory, e.g. confinement, are already contained in the limit in which the number 
of colours Nc goes large, with g'^Nc fixed. Corrections 0{1/Nc) can be treated as 
a small perturbation. A consequence of this conjecture is that also quark loops 
can be viewed as a small perturbation, apart from their effect on the scale of the 
theory. Indeed apart from the loop with two vertices which is proportional to 
g'^Nf ~ Nf/Nc ~ 1, and which enters in the (3 function, loops with more vertices 
have additional factors 1/Nc and are negligible. According to this conjecture also 
the mechanism of confinement and the corresponding order parameter should then 
be marginally affected by the presence of quarks. 



2) 



Understanding the ground state is also important to understand why perturbation 
theory works at small distances. Perturbative quantization describes interaction of 
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quarks and gluons, and the ground state is the Fock vacuum. Quarks and gluons 
are not observed in nature, and the Fock vacuum is certainly not the ground state. 
This reflects in the renormalized perturbation expansion as a lack of convergence, 
even in the sense of asymptotic expansion0. 

2. Confinement in Nature. 

Colour is confined in nature. The expected ratio of abundance of quarks Uq to abun- 
dance of nucleons Up is in the standard cosmological model 



The experimental upper limit is 

^ < 10-2^ , (2) 

Up 

coming from Millikan like experiments on ~ 1 of matter. 

The estimate (|l]) is conservative. If we assume no confinement and T is the temperature 
at which quarks decouple their effective mass is ~ T. The reactions 

g + g — s> mesons q + q ^ q + baryons 

are esothermic: let a the corresponding cross section and (Tq = lim^_»o crw. Then quarks 
will decouple when 

Uqao = G^j'^T^ . (3) 
Since ~ this implies 



G^^ 10 

— — 



1/2 Tn-18 

^ - , (4) 

(To ~ i^^"^- The ratio ([1|) corresponds to = T = lOGel^. 

The factor 10^^^ between the observation and the expectation cannot be explained 
by fine tuning of a small parameter. Like the experimental limit on the resistivity of a 
superconductor, it can only be explained in terms of symmetry. 

A suggestive idea in that direction is that vacuum is a dual superconductor p. The 
chromoelectric field between a. q q pair is constrained by dual Meissner effect into an 
Abrikosov flux tube with energy proportional to the distance. 

A relativistic version of the free energy of a superconductor, which is the analog of 
effective action in field theory, is 

G = -\f,,Fp, + {D,vr{D,v) + V{^) , (5) 
where 

Df^if = {dp, - iqA^,)if (6) 
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is the covariant derivative and 

n¥^)=/iVV-^(¥'V)' (7) 

is the effective potential, fi'^ and A are funtions of the temperature, and fi'^{T) > in the 
superconducting phase, where the potential has a mexican hat shape. 

Putting (p = /e*''^, with / > gives D^ip = — zg(y4^ — d^9)fe^'^^. Under a gauge 
transformation ^ + d^A 

d^e ^ d^e + d^K (8) 

and A^ = A^ — d^9 is gauge invariant. Moreover 

F^.u = d^A, - d,A^ = F^, (9) 
and the free energy can be rewritten as 

G = -\f,,F,, + q'fA.A, + djdj + - . (10) 
The equations of motion are 

d,F,, + q'fA, = , id' + fi')f = 2Xf (11) 



A static solution in the gauge Aq = has E = 0, f = f = ji^/2\ and eg. ([TT|) reads 

VAH + q'^pA = 0, (12) 
V^i? - q^pH = . (13) 

Eq. (|l^) means that also in the absence of electric field there is a permanent current, or 
that a = oo. Eq. ([T3|) is Meissner effect: the field A penetrates by a length (qf)^^- At 
large distance from the center of a flux tube A = or 

e^irff = e$(5) = 27m , (14) 

which is the Dirac quantization condition. Abrikosov flux tubes have monopoles at their 
ends. 

Those phenomena are a consequence of symmetry [^i the order parameter is / = (|<^|), 
or the non vanishing v.e.v. of a charged operator. For dual superconductivity the signal 
of the phase should be the v.e.v. (/x) of an operator carrying magnetic charge. 

3. Phenomenology of confinement on the lattice. 

Lattice produces evidence for confinement. Wilson loops, defined as parallel transport 
along square contours in space time, provide the static force between qq pairs, in the limit 
of large T 

W{R,T) ~ exp{-aV{R)T) (15) 

T— >oo 
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The area law observed in lattice gauge theory[0 

T) ~ exp(-ai?r) (16) 

means 

V{R) = aR (17) 

or that an infinite amount of energy is needed to pull the two particles at infinite distance 
from each other, a is the string tension, related to the slope of Regge trajectories. 

Also chromoelectric fiux tubes between qq pairs are observed, with transverse size 
~ 0.5fmp|,|^. The colour orientation of the chromoelectric field inside them can also be 
studiedi,|rg]. 

Finally the collective modes of the string formed by the fiux tube can be analysed[^]. 



All these facts support the picture of confinement as due to dual superconductivity of 
vacuum. A microscopic understanding is however needed. In particular monopoles which 
condense in the vacuum have to be identified. 

4. Monopoles 

In QED, which is a U{1) gauge theory, magnetic charges are omitted, since they are 
not observed in nature. As a consequence the general solution of Maxwell's equations can 
be given in terms of a vector potential A^. The field strength tensor 

F^u = d^A, - d,A^ (18) 

obeys the equations 

d^F^, = > . (19) 
The absence of magnetic charges indentically follows from eq. ([TE|) . The dual tensor F*^^ = 



^^^lupaF^"^ is, by virtue of eq. ([T8| ) identically conserved: 

d,F;^ = . (20) 

Eq. (|20|) is known as Bianchi identity. 

The only way to have a monopole and to preserve Bianchi identity is [|l^ to introduce 
a singularity, and consider the monopole as the end point of an infinitely thin solenoid 
(Dirac string), which can be made invisible if the parallel transport of any charge q around 
it is trivial or if 

exp (^iq (j) A - dx^ = 1 . (21) 

The line integral is intended on a path which encircles the string and is equal to the 
magnetic fiux, or to An the magnetic charge of the monopole. Eq. (^) implies 4nqM = 
2mi or qM = |, which is the celebrated Dirac quantization condition. As a consequence 
the theory becomes compact. 

In non abelian gauge theories, in the familiar multipole expansion the monopole term 
obeys abelian equation of motion, has a Dirac string and a number of independent abelian 



magnetic charges which is A^— 1 for the gauge group SU (N) . The 't Hooft-Polyakov[|l4 
[TsU monopole of the 5*0(3) Georgi-Glashow model obeys this classification. 
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5. Monopoles in QCD 

To understand the monopoles in QCD we shall phrase the classification of ref. |13[ in 
the language of ref. |T^, or in terms of "abelian projections". We shall refer to SU{2) for 
simplicity: extension to SU{N) is trivial|l7]. 

Let (f){x) ■ a be any local operator belonging to the adjoint representation. We define 
= I ^1^1 1 . (f) is well defined except at zeros of Consider the field strength 

tensor 

F^,{x) = ■ G^,{x) - ■ a Z^.0(x)) (22) 

where Gfj,i,[x) = d^Ai, — dyA^ and = — gAA is the covariant derivative. The 
coefficient of the second term in eq. (|2^) is chosen in such a way that the quadratic term 
An A Ay cancels with the first term. Both terms are gauge invariant under regular gauge 
tranformations. 

A gauge transformation U{x) which brings along the 3 axis = (0, 0, 1), and diago- 
nalizes 0(x) ■ a is called an abelian projection. After abelian projection 

F^y = d^Al - dyAl (23) 

is an abelian field. This holds in all points where U{x) is regular. Defining the dual tensor 



F^ly as 



the magnetic current j*^ is defined as 

d.K'^ = Jf (25) 

and is identically conserved. It is zero except at the singular points of U{x), where 
monopoles can appear. Thus 

= (26) 

defines a magnetic U{1) symmetry of the theory. It is not a subgroup of the gauge group 
since both F^^y and j^^^ are colour singlet. If the vacuum is not invariant under that U{1) 
symmetry, monopoles condense like the Cooper pairs and there is dual superconductivity. 
Notice that 

1. there is a magnetic U{1) symmetry for each operator 0(x) in the adjoint represen- 
tation. 

2. Under the abelian projection U{x), due to singularities, the field strength tensor 
G^y acquires a singular component 

G^y ^ UG^yU-' + . (27) 

The regular part can have monopole sources. G^*"^ describes Dirac strings starting 
from the monopoles. 
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The strategy will then be to detect condensation of different monopole species by mea- 
suring with numerical simulations across the deconfining transition a "disorder" parame- 
ter, which detects dual superconductivity. The usual order parameter is the Polyakov line. 
Our disorder parameter will be zero in the deconfined phase, where the order parameter 
is non zero, and different from zero in the confined phase, where it is zero. The concept of 
disorder parameter is typical of systems which admit a dual description They usually 
have extended structures with non trivial topology (monopoles in QCD), which condense 
in the disordered phase. In a dual description these structures are described by local 
fields, and the definition of order and disorder is interchanged. Before giving the results 
and a few details on how they have been obtained, we shall present the expectations. 



6. Expectations vs. results 

As we have seen dual superconductivity of the vacuum is not a well defined concept. 
There are infinitely many choices for the operator (j){x), and for each of them ther can or 
can not be condensation. What is the good choice, if any? 

A. There is club of practitioners of the "maximal abelian projection", saying that their 
choice is better than others. In fact with this choice the abelian field "dominates" 
the configurations, in particular the part of it which is produced by monopoles. 
This can prove convenient to attempt a construction of effective lagrangeans, but 
in principle does not preclude any pattern of symmetry. 



B. There is a conjecture that all abelian projections are equivalent P 

Most probably both attitudes reflect our imperfect knowledge of the symmetry of the 
disordered phase. 

Discriminating between (A) and (B) is possible on the lattice. The results that we have 
obtained by a systematic study of dual superconductivity show unambiguously 

1. That confinement is a transition from normal to dual superconductor ground state. 
This is strong evidence that the mechanism of confinement is indeed dual supercon- 
ductivity of the vacuum. 

2. A few different abelian projections have been analyzed. All of them show the same 
behaviour. The scenario (B) seems to be true. 

This is reassuring for the validity of the mechanism itself. If only one abelian projection 
would show dual superconductivity only the particles with non zero charge with respect 
to that U{1) could be confined: there exist states for which that charge is zero, e.g. the 
gluon which is parallel to 0. Moreover the colour direction of the electric field in the flux 
tubes observed in the lattice should also be parallel to 0, being the electric partner of the 
magnetic U{1) of the monopoles. This has been shown to be not true[^|IU|. Both these 
facts are naturally explained if the scenario (B) is at work. 

We conclude by giving a few details on the technique used to detect dual supercon- 
ductivity. The technique has been checked in many well known systems showing order 



disorder duality versus traditional descriptions[|T^]. 
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7. The disorder parameter 

An operator fi is constructed wliicli carries non zero magnetic charge with respect to 
the U{1) under study. Finite temperature is realized on the lattice by the usual ther- 
modynamical recipe of having euclidean time running from to 1/T (the temperature), 
with periodic boundary conditions for barions, antiperiodic for fermions. On a lattice this 
is done by using a size A^^^ x A^^, with A^^ > Nt and a{P)Nt = l/T {(3 = 2N/g^). By 
renormalization group arguments a(/?) ~ ^exp(— 6o/9), or T = jf^ exp{bo(3). 

The technique used to construct yU is inspired to ref. [0 and ||2^ and is a complicated 



version of the simple formula for translations in elementary quantum mechanics 

e'P-\q) = \q + a) . (28) 

In the Schrodinger representation, the field A{x,t) plays the role of x, the conjugate 
momentum n(x, t) the role of p and 



= expVJ 'j \A{x,t)) = \A{x,t) + A{x-y)) 



/i|A(f,t)) = expV ' J \A{x,t)) = \A{x,t) + A{x -y)) . (29) 

If A is the field of a monopole at y, /i is indeed the creation operator for a monopole, at 
site y and at time t. When inserted in the Feynman integral the operator is nothing but 
a linear term in the conjugated momentum added to the lagrangean and hence a shift of 
n at time t. Care is needed to adapt the definition (^) to a compact system, in a form 

which does not depend on the choice of the gauge for the classical field A. This can be 
done, and the result is, as sketched above ||T7|Jl9[ 



/ ^ Z\S + AS] , , 

with AS* different from zero on a hyperplane at constant xq. 

Being the exponential of a sum on A^^ sites (yu) is subject to strong fluctuations. Nu- 
merically it is better to measure the quantity 

P = ^ log(p) = - + A^)5+A5 (31) 



and to reconstruct (/i) as 

r/9 



(//) = exp 







p(/5') dfi' 



(32) 



A typical shape of p is shown in flg. |1|, one of (/x) in flg. 0. As is well known pT| (/x) as an 
analytic function of /3, can not vanish indentically in the deconflned phase if the number 
of degrees of freedom is flnite. An extrapolation to A^, — ^ C)0 must be done by flnite size 
analysis. Fig. |] shows that a few different abelian projections behave in the same way. 
Fig. I shows that two different monopole species of ^[/(S) also behave in the same way. 

As for the extrapolation to A^g — > oo three different ranges of /3 are considered: 
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Figure 1. p vs. P for SU{2) gauge theory. Figure 2. (/x) vs. /3 for SU{2) gauge theory. 



Plaquette projection, lattice 16 x 4. 



Plaquette projection, lattice 16 x 4. 



1. (3 ^ (3c- there perturbative theory is at work and p ^ —cNs + with c a positive 
constant. As A^s ^ oo p ^ — oo and {p) = 0. 

2. /5 -C Pc'- there p has a finite limit as Ns —>■ oo, and (p) is finite. 

3. (3 ^ Pc- there we expect (p) ~ (/?c — /?)^- From dimensional analysis 



(33) 



Where ^ > a, and if A^^ > Nt 



(p) = Ar> 0, ^, 



(34) 



and since ^ ~ (/?c — /?) we get the scaling law 



P 



(35) 



The quality of the scaling is shown in fig. ^ for SU{2) and in fig. ^ for SU{3). It 
gives a determination of the critical indices z/, 6, and of Pc- For Sf/ (2) we get 



u = 0.63(5) 

Pc = 2.30(2) AT, = 4 
6 = 0.20(8) 
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Figure 3. p vs. (3 for different abelian 
projections. SU{2) gauge theory, lattice 
12^ X 4. 



Figure 4. p vs. (3 for different abelian 
generators of SU{?>). Polyakov projection, 
lattice 12^ x 4. 



y and /3c are in agreement within errors with independent determinations, and u 
indicates a second order phase transition. 

For SUiK) 

V = 0.33(2) 

/?e = 5.70(3) Nt = A , 
S = 0.54(4) 

indicating that the transition is first order. 



The method used has been tested on a number of known systems and understood in 
its details. The result show beyond any reasonnable doubt that dual superconductivity 
occurs, in different abelian projections, in connection with confinement. 

The part of the work reported due to our group has been done in collaboration with 
L. Del Debbio, G. Paffuti, P. Fieri, B. Lucini, D. Martelli in the last few years. Their 
contribution was determinant to the results. 
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